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Abstract 

We consider the production of IU7 and Z 7 pairs at hadron colliders. We report on 
the complete fully differential computation of radiative corrections at next-to-next- 
to-leading order (NNLO) in QCD perturbation theory. The calculation includes the 
leptonic decay of the vector boson with the corresponding spin correlations, off shell 
effects and final-state photon radiation. We present numerical results for pp collisions 
at 7 and 8 TeV and we compare them with available LHC data. In the case of Z7 
production, the impact of NNLO corrections is generally moderate, ranging from 8% to 
18%, depending on the applied cuts. In the case of IU7 production, the NNLO effects 
are more important, and range from 19% to 26%, thereby improving the agreement of 
the theoretical predictions with the data. As expected, the impact of QCD radiative 
corrections is significantly reduced when a jet veto is applied. 
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1 Introduction 


The discovery of a new scalar resonance in the search for the Standard Model (SM) Higgs boson [1] 
2] is a milestone in the LHC physics programme. The properties of this new particle closely 
resemble those of the Higgs boson, but further work is needed to clarify if it is really the Higgs 
boson predicted by the SM, or something (slightly) different. Vector-boson pair production has a 
prominent role in this context. It represents an irreducible background to Higgs and new-physics 
searches, and, at the same time, it provides information on the form and the strength of the 
vector-boson gauge couplings. The interactions of W and Z bosons with photons are particularly 
interesting as they test the WW^ and ZZ'y couplings, which are predicted by the non-Abelian 
SU( 2 )l <g) U(1) Y gauge group. 

Constraints on WW^ and ZZ 7 anomalous couplings have been obtained at LEP [3]. At 
hadron colliders, studies of V 7 final states have been first carried out at the Tevatron M, and 
they were used to set more stringent limits on anomalous couplings. The high-energy proton- 
proton collisions at the LHC allow us to explore the production of V 7 ( V = W ± ,Z) pairs in a 
new energy domain. Measurements of V 7 final states have been carried out by ATLAS [Mil] and 
CMS mm using the data sets at centre-of-mass energy = 7 and 8 TeV. These measurements 
have been compared to the SM predictions and used to improve the limits on anomalous couplings 
and on the production of possible new resonances. 

When considering the V 7 final state, besides the direct production in the hard subprocess, the 
photon can also be produced through the fragmentation of a QCD parton, and the evaluation of 
the ensuing contribution to the cross section requires the knowledge of a non-perturbative photon 
fragmentation function, which typically has large uncertainties. The fragmentation contribution 
is significantly suppressed by the photon isolation criteria that are necessarily applied in hadron- 
collider experiments in order to suppress the large backgrounds. The standard cone isolation, 
which is the standard choice in the experiments, suppresses a large fraction of the fragmentation 
component. The smooth cone isolation completely suppresses the fragmentation contribution [IB] , 
but the algorithm is difficult to be implemented experimentally. 

The present status of theoretical predictions for V 7 production at hadron colliders is as fol¬ 
lows. The V'f cross section is known in next-lo-leading-order (NLO) QCD [T6llT7j . and the leptonic 
decay of the vector boson has been included in Ref. [T 8 ] . In the case of Z 7 the loop-induced gluon 
fusion contribution, which is formally next-to-next-to-leading order (NNLO), has been computed 
in Ref. [13,120], and the leptonic decay of the Z boson, together with the gluon-induced tree-level 
NNLO contributions, have been added in Ref. [2TJ . The NLO calculation for V 7 , including photon 
radiation from the final-state leptons, the loop-induced gluon contribution and the photon frag¬ 
mentation at LO have been implemented into the general purpose numerical program MCFM [22]. 
Electroweak (EW) corrections to V 7 production have been computed in Refs. [23121] • The full 
NLO EW corrections to IL 7 production with leptonically decaying W bosons, taking into ac¬ 
count all off-shell effects in the complex-mass scheme, and all effects originating from initial-state 
photons, have been computed in Ref. |25j . For W 7 production, the NLO computation has been 
matched to a parton shower according to the MiNLO prescription [26] in Ref. m 

In a previous letter [28] we have presented the results of the first complete NNLO calculation 
for Z 7 production. In the present paper, we extend this calculation to the complete class of V 7 
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production with leptonic decays, namely to both Z 7 production with visible (Z —> £ + £~) and 
invisible (Z —> vv) Z- boson decays, and to IF 7 production with the respective decays W + —y u£ + 
and bb - —> £~u. Off-shell effects and final-state photon radiation are consistently includecfj]. For 
these production channels, we present detailed results on fiducial cross sections and distributions 
at yfs = 7 and 8 TeV, and provide comparisons to ATLAS data, where available. 

The paper is organized as follows. In Section [2] we provide some technical details of our compu¬ 
tation and discuss the particular challenges in the cancellation of infrared singularities. Section [3] 
contains our theoretical predictions for all V 7 processes as well as a comparison with experimental 
data. We consider Z 7 production in the visible (Z —y £ + £~) and invisible (Z —y vv) decay chan¬ 
nels in Section [3731 and Section [TT and W 7 production in Section [3751 In Section [3761 we discuss 
the different impact of QCD radiative corrections in the W 7 and Z 7 processes and its physical 
origin. In Section [4] we summarize our results and comment on the remaining uncertainties. 


2 Details of the calculation 


In this Section we discuss the details of our calculation. We first point out that the notation 
“V 7 ” suggests the production of an on-shell vector boson plus a photon, followed by a factorized 
decay of the vector boson. Instead, we actually compute the NNLO corrections to the processes 
pp —* £ + £~^f + X, pp —* vivff + X, and pp —* £v^ + X, where, in the first case, the lepton pair 
£ + £~ is produced either by a Z boson or a virtual photon. All contributions where the final-state 
photon is radiated off the charged leptons are consistently included (see Figure [Hand Figure [2]). 
The shortcuts “Z 7 ” and “IF 7 ” are used only for convenience. 



(a) Topology la (b) Topology lb (c) Topology Ha 

Figure 1: Feynman diagrams contributing to Z 7 production at 


(d) Topology lib 
Born level. 



(a) Topology la (b) Topology lb (c) Topology Ha (d) Topology III 

Figure 2: Feynman diagrams contributing to IF 7 production at Born level. 

The NNLO computation requires the evaluation of tree-level scattering amplitudes with up to 
two additional (unresolved) partons, of one-loop amplitudes with up to one additional (unresolved) 

tFirst results from this calculation on thy production have been presented in Ref. [29) . 
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parton [5011ST] , and of one-loop squared and two-loop corrections to the Born subprocess (qq —> 
£ + £ - 7 and qq —> ugi?g 7 for Z 7 , qq' —> fo'g'y for W 7 ). Furthermore, processes with charge-neutral 
final states receive loop-induced contributions from the gluon fusion channel (gg —> £ + £~ 7 and 
gg —* viva). In our computation, all required tree-level and one-loop amplitudes are obtained from 
the OpenLoops generator [32^ . which implements a fast numerical recursion for the calculation 
of NLO scattering amplitudes within the SM. For the numerically stable evaluation of tensor 
integrals we rely on the COLLIER library [33], which is based on the Denner-Dittmaier reduction 
techniques |34tl35] and the scalar integrals of [36]. 

The two-loop corrections to the Drell-Yan-like Born processes, where the photon is radiated 
off the final-state leptons, have been available for a long time [37]. The last missing ingredient, 
the genuine two-loop corrections to the V 7 amplitudes, have been presented in Ref. [38]. 

In Section 12.11 we give a sketch of the q y subtraction formalism, the method we use to combine 
the (separately divergent) ingredients of the NNLO calculation to obtain quantitative predictions. 
The bookkeeping of all partonic subprocesses and the numerical integration of the different cross 
section contributions is managed by the fully automatized Munich framework [§, which is described 
in Section 12.21 In Section 12.31 the double-virtual and gluon fusion contributions are discussed, and 
Section I2A\ addresses the real-emission part of the NNLO cross section with a discussion of the 
numerical stability issues, namely the reliable evaluation of one-loop amplitudes in the real-virtual 
part and the numerically stable phase-space integration in the double-real part. 


2.1 The qt subtraction formalism 


The implementation of the various scattering amplitudes in a complete NNLO calculation is a 
highly non-trivial task due to the presence of infrared (IR) singularities at intermediate stages of 
the calculation, which prevents a straightforward implementation of numerical techniques. Various 
methods have been proposed and used to overcome this difficulty. They are based either on 
extensions of the subtraction method [391142] at NNLO [431146] . or on sector decomposition M- 
More recently, also a combination of the subtraction method with sector decomposition has been 
proposed [491150] . The q t subtraction formalism [5T] is an independent method to handle and cancel 
IR singularities at the NNLO. In its present formulation the method applies to the production 
of a colourless high-mass system F in generic hadron collisions, and has been applied to the 
computation of NNLO corrections to several hadronic processes PBIEHM]. According to the q y 
subtraction method [51], the pp —> F + X cross section at (N)NLO can be written as 


dcr 


F 

(N)NLO 


njF 

rt (N)NLO 


® dcr LO 


+ 


d(T (N)LO - d(T 


CT 

(N)NLO > 


( 1 ) 


where dcr^^o represents the cross section for the production of the system F plus one jet at 
(N)LO accuracy, and can be evaluated with any available NLO subtraction formalism. The (IR 
subtraction) counterterm dcr^ NLO is obtained from the resummation of logarithmically-enhanced 
contributions to gy distributions [58]. The practical implementation of the contributions in the 
square bracket in Eq. (JT|) is described in more detail in Section 12.41 

*The OpenLoops one-loop generator by F. Cascioli, J. Lindert, P. Maierhofer and S. Pozzorini is publicly 
available at http: //openloops . hepf orge. org 

§ Munich is the abbreviation of “MUlti-chaNnel Integrator at Swiss (CH) precision”—an automated parton 
level NLO generator by S. Kallweit. In preparation. 
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The ‘coefficient’ ”H(n) NLO enc °des the loop corrections to the Born-level process and also com¬ 
pensated for the subtraction of dcr^ NL0 . It is obtained from the (N)NLO truncation of the 
pro cess-dependent perturbative function 

M p = l + ^M F(I >+f^')V (2) + ... . (2) 

7T V 7T / 

The NLO calculation of da 1, requires the knowledge of T-L F ^\ and the NNLO calculation also 
requires The general structure of H b(1 ^ is known [59]: 'H hFL ' 1 is obtained from the process- 

dependent scattering amplitudes by using a process-independent relation. Exploiting the explicit 
results of for Higgs [60] and vector boson [61] production, the process-independent relation 
of Ref. [59] has been extended to the calculation of the NNLO coefficient [62]. These 

results have been confirmed through an independent calculation in the framework of Soft-Collinear 
Effective Theory (SCET) [63] [64]. The counterterm dcr^ NLO only depends on i.e. for a 

NNLO computation, it requires only as input, which can be derived from the one-loop 

amplitudes to the involved Born subprocesses. 

As can be seen from the discussion above, a significant part of the ingredients needed to perform 
the NNLO computation are actually NLO-like in nature, allowing us to build the implementation 
of our NNLO calculation upon existing NLO tools. We have based this calculation on the Munich 
framework, which is briefly discussed in the following. Necessary extensions of this framework are 
addressed as well in the rest of this Section. 


2.2 Organization of the calculation within the Munich framework 

MUNICH is a fully automatized framework for the computation of fixed-order cross sections for 
arbitrary SM processes up to NLO accuracy, written in C++. After the hadronic process has 
been specified, Munich takes care of the bookkeeping of all partonic subprocesses that need 
to be taken into account. It automatically generates adequate phase-space parametrizations for 
each partonic subprocess by exploiting the resonance structure of the underlying (squared) tree- 
level Feynman diagrams. These parametrizations are combined using a multi-channel approach 
to simultaneously flatten the, in general complicated, resonance structure of the amplitudes and 
thus guarantee a reasonable convergence of the numerical integration. Several improvements like 
an adaptive weight-optimization procedure are implemented as well. 

Munich was originally developed for the NLO calculations of [65bl67] . where only massless 
colour-charged particles were involved. To account for the mediation of IR singularities between 
the different phase spaces of virtual and real corrections, the Catani-Seymour dipole subtraction 
method [TT[ I42] was implemented. In Ref. [68], the framework was extended to massive QCD 
according to Ref. [69], and in Ref. m the extension to complete SM calculations, i.e. including 
also EW corrections at NLO, was presented. All dipole terms necessary to numerically cancel 
the soft and collinear divergences of the real corrections are constructed at runtime from spin- 
and colour-correlated matrix elements of the underlying partonic Born subprocesses. Moreover, 

^More precisely, while the behavior of dcr^ NLO as —> 0 is dictated by the singular structure of , 

its non-divergent part in the same limit is to some extent arbitrary, and its choice determines the explicit form of 

n_/F 

n '(N)NLO‘ 
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additional phase-space parametrizations based on the reduced dipole kinematics are generated 
and included in the multi-channel, supplementing the parametrizations based on real-emission 
kinematics. Analogously, the analytically integrated dipoles, the so-called K + P terms and the 
/-operator, which compensate for the dipole terms subtracted on the real-emission side, are auto¬ 
matically constructed at runtime from colour-correlated matrix elements of the Born subprocesses. 
For the evaluation of all involved matrix elements up to the one-loop order, MUNICH provides an 
automatic interface to amplitudes generated by the OpenLoops generator [32] . 

The guiding principle in our NNLO implementation is to keep the additionally introduced pro¬ 
cess dependence to the bare minimum, essentially limiting it to the two-loop amplitudes entering 
7-L h in Eq. ([I]). All other process dependent information entering the various pieces in Eq. ([Tj) 
have been expressed in terms of NLO quantities available inside Munich via its interface to 
OpenLoops. 


2.3 Double-virtual and gluon fusion contribution 

In our implementation of the q t subtraction formalism, the only component that needs to be 
provided on a process-by-process basis are the two-loop—and one-loop-squarecH— amplitudes to 
the Born processes, which enter the coefficient , H( ? N ) NLO in Eq. (JTJ). Ref. [3S] provides the analytical 
expressions of the resonant hclicity amplitudes for Eq production (e.g. qq —?• Zq — > / + / - q) in 
terms of two-dimensional harmonic polylogarithms up to weight 4. The non-resonant (final-state 
radiation) contribution (e.g. qq —* £ + £~ —> / + /“q) is described by the two-loop quark form factor 
from Ref. [37], We have implemented the hclicity amplitudes directly into a C++ library. For 
the numerical evaluation of the harmonic polylogarithms we use the tdhpl library mi Our 
implementation allows for several thousand amplitude evaluations per minute, rendering the time 
spent on computing the two-loop contribution negligible compared to the time needed for the real 
emission contribution. 

For the processes with an electrically neutral final state at Born level, i.e. pp —> / + /“q and 
pp —* z+Xqq, the loop-induced gluon fusion processes gg —$■ Z + Z“q and gg —$■ iqX+q represent 
a separately IR-finite and gauge-invariant part of the cross section. Though of O (a§), and thus 
formally of NNLO, the gluon fusion contribution is often included already at NLO (see for example 
MCFM [22].) The reason for that is the assumption that the large gluon luminosities at the LHC 
could compensate for the additional as suppression, and hence the contribution could numerically 
become as important as the NLO corrections themselves. As we will see in Section 13.31 and 
Section l3~4l this is not the case in Zq production. The gluon fusion contribution is small even 
when compared to the other NNLO corrections. In our computation, these (finite) one-loop- 
squared gluon fusion amplitudes are evaluated using OpenLoops. 

II The OpenLoops generator also provides the non-finite one-loop-squared matrix elements in terms of the 
corresponding coefficients of the Laurent series. These results are, however, only used as a numerical check in this 
calculation. 
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2.4 Real-emission and counterterm contribution 


All NNLO contributions with vanishing total transverse momentum q x of the final state system 
F are collected in the coefficient which is discussed in the previous Section. The remain¬ 

ing part of the NNLO cross section, namely the difference in the square bracket in Eq. (JTi) , is 
formally finite in the limit gx —t 0, but the terms separately exhibit logarithmic divergences in 
this limit. The counterterm dcr^ NLO is integrated over the n-particle Born phase-space**!. while 

the term da^Lo* i n Eq. ([I]), called the real-emission contribution in the following, involves two 
different phase spaces with one and two additional QCD partons, respectively, and thus needs 
to be integrated separately over the respective (n + 1)- and (n + 2)-particle phase spaces. To 
achieve a numerical cancellation of the singularity, a technical cut on q £ needs to be introduced to 
render both terms separately finite. In this sense, the q t subtraction method works very similarly 
to a phase-space slicing method at NLO. In practice, a technical cut, r cut , on the dimensionless 
quantity r = q^/M, where M denotes the invariant mass of the final-state system F, turns out to 
be a more convenient choice than a cut on q x itself. 

By construction, the counterterm dcr^ NLO cancels all gx-divergent (logarithmic) terms from 
the real-emission contributions, implying that the r cut dependence of their difference does not only 
vanish in the limit r cut —)■ 0, but should become arbitrarily small for sufficiently small values of r cut . 
In practice, however, as both the counterterm and the real-emission contribution grow arbitrarily 
large for r cut —>■ 0, the statistical accuracy of the Monte Carlo integration degrades, preventing one 
from pushing r cut arbitrarily low. In general, the absence of any strong residual r cut dependence 
in the difference between the real contribution da^^ and the counterterm dcr^ NLO provides 
a strong check on the correctness of the computation since any (significant) mismatch between 
the contributions would result in a divergent cross section in the limit r cut —>■ 0. To monitor the 
r cut dependence without the need of several CPU-intensive runs, our implementation allows for 
simultaneous cross section evaluations at arbitrary r cut values. 

Typical values of r cut are of the order of 10 -3 , requiring a numerically stable evaluation of 
dffivLO^ down to the IR-divergent region where the transverse-momentum of the Vj system 
reaches values of about 0.1 GeV. While the IR-enhanced phase-phase region with a very-low-gx 
jet is excluded in standard NLO calculations, it needs to be resolved at very high precision in a 
NNLO computation. This issue obviously challenges a dedicated NLO program like Munich. We 
have slightly modified the phase-space generation to achieve the required precision and reliability 
of our NNLO results within the MUNICH framework. 

The real-emission contribution to the NNLO cross section is, apart from the presence of a 
very-low-g T jet discussed above, a usual NLO calculation, and can thus be treated with the 
Catani-Seymour dipole subtraction method implemented within MUNICH. Consequently, it can 
be split into a real-virtual (RV), a real-collincar (RC) and a double-real (RR) contribution, 

d(J NLO Ct = dff RV + d(J RC + d(J RR- ( 3 ) 

Here, the mediation of all IR divergences (given a finite r cut value) by dipole subtraction is im¬ 
plicitly understood, i.e. each contribution on the right-hand side is finite and can be numerically 
integrated over the respective phase space. 

**The integration over qx only acts on an explicit <jx dependence of the phase-space weight. 
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Whereas the real-collinear sub contribution does not exhibit any peculiar issues, the other two 
subcontributions on the right-hand side pose different challenges when integrated into the deep 
IR (< 7 t > 0) region. The real-virtual contribution requires the evaluation of one-loop matrix 
elements to Vy + jet production, which has become a standard task for automatic one-loop tools. 
However, as we evaluate the matrix elements in the deep IR region, i.e. far away from the phase- 
space region that is relevant for 2 —> 3 hard scattering processes at NLO, numerical instabilities 
in the amplitude evaluation might be a source of concern. To address this issue, OpenLoops 
implements a fully automated system that monitors the numerical accuracy of loop amplitudes 
and cures possible instabilities at runtime. This system exploits the fact that, for the reduction 
to scalar integrals, OpenLoops allows one to flexibly switch from tensor-reduction [[3TI135] to OPP 
reduction m algorithms. 

To perform the reduction from tensor to scalar integrals, we make use of the Collier li¬ 
brary [33] that consists of two independent implementations of the Denner-Dittmaier reduction 
algorithm [34l[35]- The presence of potential residual instabilities is tested by comparing the two 
implementations of tensor reduction. In presence of instabilities that exceed one permille of the 
Born amplitude, the one-loop amplitude is reevaluated using the CutTools m implementation 
of OPP reduction in quadruple precision. In this case OneLOop m is used for the evaluation 
of scalar integrals. Finally, the accuracy of the result in quadruple precision is verified by a con¬ 
sistency check based on the rescaling of all dimensionful variables. For all processes considered 
in this paper we find that the Denner-Dittmaier reduction works fast and reliably for more than 
99% of the phase space points. This allows one to restrict the usage of quadruple precision to a 
tiny fraction of points with a minor impact on the total runtime of the calculation. 

The double-real contribution, on the other hand, only involves tree-level amplitudes, but con¬ 
tains an additional unresolved parton in the final state. Moreover, it involves several dipole terms 
located on different (n + l)-particle phace spaces. Like any other phase-space cut, the restriction 
r > r cut has to be applied on the respective phase space and can thus lead to miscancellation 
issues, which are well-known from NLO calculations, but might be amplified here when happening 
in the deep IR region. These complications render the numerical phase-space integration in the 
deep IR region more challenging compared to the integration of the single-emission phase space. 
Munich already provides a state-of-the-art multi-channel phase-space integrator, which greatly 
speeds up computations compared to a classical VEGAS integration procedure, and guarantees a 
reasonably stable and reliable convergence behaviour. In the deep IR regions, however, even an 
advanced multi-channel approach might fail to capture all relevant contributions. We have thus 
extended the Munich integrator by an additional VEGAS-like importance sampling on top of the 
multi-channel parametrization. This hybrid approach results in an extremely efficient and reliable 
phase-space integration. 

The numerical information on the r cut dependence of the cross section is used to perform an 
extrapolation to r cut = 0, which can in turn be used to quantify the uncertainty due to the hnite 
r cut value. This estimated uncertainty is combined with the usual statistical integration error to 
provide an overall estimate of the numerical precision of our NNLO prediction. For the processes 
considered in this paper, the extrapolation at r cut —> 0 turns out to be non-trivial, due to the 
interplay with the photon isolation. Nonetheless, the procedure allows us to control integrated 
NNLO cross sections to better than 1%. 


7 


3 


Results 


3.1 Setup 


For the electroweak couplings we use the so-called G M scheme, where the input parameters are 
Gf, mw, m z- hi particular we use the values Gf = 1.16639 x 10 -5 GeV“ 2 , mw = 80.399 GeV, 
mz = 91.1876 GeV, Yz = 2.4952 GeV and Y\y — 2.1054 GeV. We set the CKM matrix to unity. 
We use the MMHT 2014 [75] sets of parton distribution functions (PDFs), with densities and as 
evaluated at each corresponding order (i.e., we use (n + l)-loop as at N n LO, with n — 0,1,2), 
and we consider Nf = 5 massless quarks/antiquarks and gluons in the initial state. 

The default renormalization (p R ) and factorization (p F ) scales are set to p R = p R = Pq = 
V m v + (Pf) 2 - An estimate of missing higher-order contributions is obtained by performing p R 
and p R variations by a factor of two around the central value. We find substantial cancellations 
between renormalization and factorization scale variations in some of the calculations we are going 
to present if the restriction p R = p F is imposed. These cancellations are assumed to be purely 
accidental. To accomodate this well-known feature, we consider also antipodal variations of the 
two scales [22], i.e. setting p R = £/z 0 , Pf = Po/£ and varying £ between | and 2. In summary, we 
estimate scale uncertainties by varying p F and p R simultaneously and independently in the range 
0.5/i 0 and 2 p 0 with no constraint on the ratio Pf/Pr- 


The present formulation of the qp subtraction formalism [51] is limited to the production 
of colourless systems F and, hence, it does not allow us to deal with the parton fragmentation 
subprocesses. Therefore, we consider only direct photons, and we rely on the smooth cone isolation 
criterion [15]. Considering a cone of radius r = \J (A rf) 2 + (A0 ) 2 around the photon, we require 
that the total amount of hadronic (partonic) transverse energy E T inside the cone is smaller than 
E“ ax (r), 


£” ax (r 


7 

€ 7 pJ^ 


1 — cos r\ n 
1 — cos R) 


(4) 


where is the photon transverse momentum; the isolation criterion Er < Fi™ ax (r) has to be 
fulfilled for all cones with r < R. All results presented in this paper are obtained with e 7 = 0.5, 
n = 1 and R = 0.4. 


3.2 Comparison to experimental data 

The smooth cone isolation prescription adopted in our calculation is not yet implemented in exper¬ 
imental analyses. Measurements are usually performed by using a fixed cone isolation prescription 
(which corresponds to Eq. (J4J) with n — 0), and thus our isolation prescription does not exactly 
correspond to what is done in the experiment. However, the parameters e 7 and R needed to spec¬ 
ify the smooth cone have natural counterparts in the definition of the fixed cone, while the precise 
choice of the smoothing function (in our case parametrized by n) does only have a mild impact on 
the final result. Furthermore, recent studies carried out in diphoton production [76] show that for 
sufficiently tight isolation parameters, smooth and hard cone isolation yield very similar results. 
For the processes in this paper, we verified at NLO that the difference between using smooth and 





hard cone isolation is at the 1 — 2% levefd], i.e. well below the current experimental uncertainties 
and still smaller than the remaining theoretical uncertainties. We can thus safely compare our 
theoretical predictions with experimental data. 

Since the first results of our work have appeared [281(29], we have provided numerical predic¬ 
tions for Z 7 production to the CMS collaboration ra, and for W 7 production to the ATLAS 
collaboration m, as a background in the H —>• WW analysis. These predictions were obtained 
by using the experimental cuts adopted in the corresponding analyses. In the present paper, we 
limit ourselves to compare our predictions to the ATLAS results for W 7 and Z 7 at 7 TeV 0. 

3.3 pp —> 

In our calculation of pp —$■ £ + £~^f at y/s = 7 and 8 TeV we adopt the selection cuts used by the 
ATLAS collaboration [9], summarized in Tabled! We require the photon to have a transverse 
momentum of > 15 GeV (soft p^ cut) or p\ > 40 GeV (hard cut) and pseudorapidity 
| 7 7 | < 2.37. Each of the charged leptons is required to have p ^ > 25 GeV and \rj e \ < 2.47, and 
the invariant mass of the lepton pair must fulfill mg+g- > 40 GeV. We require the separation in 
rapidity and azimuth A R between the leptons and the photon to be AR(£, 7 ) > 0.7. Jets are 
reconstructed with the anti- fop algorithm [78] with radius parameter D = 0.4. A jet must have 
p > 30 GeV and | 7 jet | < 4.4. We require the separation A R between the leptons (photon) and 
the jets to be AR(£/'j,]et) > 0.3. At yfs = 8 TeV, the jet definition is slightly adjusted by using 
| 7 jet | < 4.5 instead of \rf et \ < 4.4, adapting to the ATLAS Run II standard. With respect to 
resolved jets in the final states, we will consider both the inclusive (Vj et > 0) and the exclusive 
(A jet = 0) case. 

The predicted cross sections with the soft p ^ cut, including the theoretical uncertainties from 
scale variations obtained as described at the beginning of Section l3Tl can be found in Table [2] 
When going from NLO to NNLO the cross section increases by 8 % (3%) in the inclusive (exclusive) 
case, respectively. The fiducial cross sections measured by ATLAS at 7 TeV [9j are also reported 
in Table (21 Both the NLO and NNLO predictions are in agreement with the experimental result, 
and the NNLO corrections improve the agreement, especially in the inclusive case. 

The reduced impact of QCD radiative corrections when going from the inclusive (A’j et > 0) to 
the exclusive (Aj et = 0) case is a well known feature in perturbative QCD calculations [79] . A 
stringent veto on the radiation recoiling against the Z 7 system tends to unbalance the cancella¬ 
tion between positive real and negative virtual contributions, possibly leading to large logarithmic 
terms. The resummation of these logarithmic contributions has been the subject of intense theo¬ 
retical studies [80b82] . especially in the important case of Higgs boson production. The reduced 
impact of radiative effects in the presence of a jet veto is often accompanied by a reduction of scale 
uncertainties. In the present case, since we are considering a process initiated by quark-antiquark 
scattering, the impact of radiative corrections is smaller than in Higgs boson production. How¬ 
ever, a reduction of scale uncertainties from the Aj et > 0 to the Aj et = 0 case is already visible 
in Table (21 and may signal the need of more sophisticated (conservative) methods to estimate 
perturbative uncertainties [801183]. 

^Obviously, the agreement also significantly depends on the fragmentation function used when employing the 
hard cone isolation, which typically has large uncertainties. 
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y/s = 7 Te V 


yfs = 8 TeV 


Leptons 

p f T > 25 GeV 
\rj\ < 2.47 

Photon 

p^ > 15 GeV (soft Pt cut) or p^ > 40 GeV (hard p ^ cut) 

|p 7 | < 2.37 

Frixione isolation with e 1 = 0.5, R = 0.4, n = 1 

Jets 

anti-/cx algorithm with D = 0.4 
pf > 30 GeV 

rf et < 4.4 |r/ et | < 4.5 

iVjet > 0 (inclusive) or ATj et = 0 (exclusive) 

Separation 

rri£+i- > 40 GeV 

AR(£, 7 ) > 0.7 

Ai?(£/ 7 , jet) > 0.3 


Table 1: Event selection criteria used in the pp £ + £ 7 analysis. 


^ [TeV] 


CAo [pb] 

C^NLO [Pb] 

0NNLO [pb] 

ATLAS [pb] 

7 

Aj e t A 0 

Aj e t = 0 

0-814915S 

1 000+4.2% 

i.zzz_ 53% 

1-0313S 

1-320^:11 

1.059+S 

±0.02 (stat) 

1.31 ±0.11 (syst) 

±0.05 (lumi) 
±0.02 (stat) 

1.05 ±0.10 (syst) 

±0.04 (lumi) 

8 

Ajet A 0 
Ajet = 0 

0.92447s 

1.3877;“ 

1 1 c'7+2.6% 
L ‘ 10 -4.5% 

1 504 +L3% 

1 1 oq+0-8% 



Table 2: pp —y £ + £~7 cross sections with the soft p ^ cut (p^ > 15GeV). Scale uncertainties are 
obtained from independent variations of pn and p,p around the central scale /x 0 , as described in 
Section [BTT1 The numerical uncertainty of the NNLO prediction from statistical error and finite 
r cut are conservatively estimated to be about 0.3%. The last column provides the corresponding 
results by ATLAS. 


^ [TeV] 


CAo [fb] 

cnlo [fb] 

onnlo [fb] 

7 

CD 

IV 

0 

73 fil+ 3 - 4% 
/o.oi_ 45% 

132 0 +4 - 2% 
ioz.u_ 40% 

154.3+ 3 ;S 

8 

J<S 

CD 

IV 

0 

84 f)Q +4 - 3% 

°^- Uy -5.5% 

1 CO 1 ±4.6% 
100 - 1 -4.5% 



Table 3: pp —> £ + £~^ cross sections with the hard p\ cut (p\ > 40GeV). Scale uncertainties are 
computed as in Table [2] The numerical uncertainty of the NNLO prediction from statistical error 
and finite r cut is conservatively estimated to be about 0 . 6 %. 
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Figure 3: Photon transverse momentum distribution in the inclusive (left) and exclusive case 
(right) at NLO (red, dashed) and NNLO (green, solid) compared to ATLAS data. The lower panel 
shows the data/theory ratio for both theory preditions. In the upper panel, only experimental 
uncertainties are shown. The lower panel also shows theoretical uncertainty estimates from scale 
variations. 


Beyond the cross section in the fiducial region, ATLAS has also provided the measured cross 
section differential in the photon transverse momentum. A comparison of the resulting distribution 
with our theoretical NLO and NNLO predictions is displayed in Figure [H] for both the inclusive and 
the exclusive case. In particular at transverse momenta pj. < 100 GeV, the inclusion of NNLO 
corrections tends to improve the agreement between data and theory. The comparison of the 
theoretical predictions to the data in the high transverse-momentum region > 100 GeV is quite 
delicate. First, the experimental uncertainty in this region is quite large. Then, EW corrections are 
expected to become sizable and negative due to large Sudakov logarithmic contributions [25|E3], 

In Figure Owe compare the NLO and NNLO predictions for the invariant-mass distribution 
of the £ + £~7 system with the distribution provided by ATLAS in Ref. [9]. For this measurement, 
ATLAS increases the transverse momentum cut on the photon to p\ > 40 GeV: the corresponding 
cross sections are reported in Table [3] The relative impact of radiative corrections is 79% and 
17% when going from LO to NLO and from NLO to NNLO, respectively. We conclude that the 
corrections are significantly larger compared to the case in which the soft p^, cut (p^, > 15 GeV) 
is applied. As the differential cross section in Figure [4] is normalized by the fiducial cross 

section, sizeable NNLO corrections are visible only in the first bin, where the agreement with 
data is slightly improved. This implies that the NNLO/NLO ratio is almost constant for larger 
invariant masses. 

The more pronounced higher-order corrections in the case in which a hard p vj) cut is applied 
can be understood by studying the £ + £~-y invariant mass distribution in a finer binning, which is 
shown in Figure [5] for both the soft and the hard p\ cuts. When the soft p^ cut is applied, the 
relative impact of the NNLO corrections is small in the region around the Z —>■ £ + £~7 peak, where 
the fiducial cross section receives its dominant contribution, and then slowly increases with the 
invariant mass. When the hard p^ cut is applied, the Z —> £ + £~7 peak is not populated at all at 
LO as the applied cuts produce a lower bound at m ^+£- 7 « 97 GeV in LO kinematics. The region 
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Figure 4: The invariant-mass distribution of the £ + £~^ system at NLO (red, dashed) and NNLO 
(green, solid), normalized to the fiducial cross section calculated at the repsective order, is com¬ 
pared to ATLAS data. The lower panel shows the data/theory ratio. Only experimental uncer¬ 
tainties are shown. 


below the boundary contributes sizably to the cross section, but is only populated beyond LO, 
i.e. in this region the NLO computation provides actually only a LO prediction. Hence the NNLO 
predictions effectively correspond to the first perturbative correction, with a comparably large K 
factor of about 1.4. The lower bound on m £+£- 7 for LO kinematics also exists with the soft p 7 
cut, namely at m £+^- 7 ~ 66 GeV. However, in this case the Z —> £ + £~^ peak is populated already 
at LO, and the region below the cut does not significantly affect the fiducial cross section. As 
already expected from Figure [H the NNLO/NLO ratio in the hard p\ case is almost independent 
of m £+^- 7 above m ^- 7 ~ 140 GeV. 

Figure O also shows the contribution from the loop-induced gluon fusion process, which, as 
explained in Section [21 respresents a finite and gauge invariant subcontribution to the full NNLO 
result. This contribution is often argued to be potentially sizable due to the large gluon luminosi¬ 
ties at the LHC. In our calculation, however, the gluon fusion contribution turns out to be small. 
It amounts only to around 6(9)% of the full O («g) correction and, correspondingly, to less than 
1 ( 2 )% of the total fiducial cross section in the soft and the hard pp case, respectively. 

3.4 pp -A 

In the pp —y Zy —y zyTyy analysis for proton-proton collisions at yfs = 7TeV, we again use the 
selection criteria applied by ATLAS [9J: compared to the pp —$■ £ + £~^ analysis, the transverse- 
momentum cut on the photon is made harder (p^> >100 GeV), and a cut on the missing transverse 
momentum, i.e. the vectorial sum of the neutrino momenta, p™ > 90 GeV, is imposed. The 
jet algorithm, the photon isolation and all other event-selection criteria are the same as in the 
pp —> £ + £~7 analysis, if applicable. 

In the yfs = 8 TeV analysis, both the photon transverse-momentum and the missing transverse- 
momentum cuts are increased to pp > 130 GeV and p™ > 100 GeV, respectively. As in the 
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Figure 5: Invariant mass distribution of the l + t~ r y system at LO (blue, dotted), NLO (red, dashed) 
and NNLO (green, solid) for the setup with > 15 GeV (left) and the setup with p\ > 40 GeV 
(right). The loop-induced gluon fusion contribution is also shown (pink, dash-dotted). The lower 
panel shows the NNLO/NLO ratio. 


Z'y —y £ + £ 7 analysis the rapidity acceptance for jets is slightly increased to | 7 /j et | < 4.5. The cuts 
are summarized in Table [4] 

The predicted cross sections at LO, NLO and NNLO can be found in Table Ej] The results 
are presented for a single neutrino species and thus have to be multiplied by a factor of three to 
obtain the complete uu'-y cross section. In the inclusive case, i.e. for iVj et > 0, we find relatively 
large NLO corrections of around 57% and 68 % and NNLO corrections of around 12% and 14% 
at y/s = 7TeV and y/s = 8 TeV, respectively. The inclusive NNLO cross section prediction at 
y/s = 7TeV is in good agreement with the cross section measured by ATLAS. In the exclusive 
case, IV^ = 0, the NNLO corrections are very small, and the scale uncertainties are reduced 
down to the 1% level. We observe quite a significant discrepancy with respect to the ATLAS 
measurement for y/s = 7TeV. The origin of this discrepancy is unclear to this point. As mentioned 
in Section 13.31 the stability of the fixed order calculation when a jet veto is applied is challenged 
and the perturbative uncertainties we find through scale variations are likely to be underestimated. 
We also point out that our NLO prediction differs significantly from the NLO prediction reported 
in Table VII of Ref. [9], that, even using MCFM, we are not able to reproduce. This could be 
due to a large parton-to-particle correction applied on the NLO result quoted by ATLAS. 

Figure [ 6 ] shows the photon transverse-momentum and the missing transverse-momentum dis¬ 
tributions. These distributions are identical for Born kinematics due to momentum conservation, 
so the difference results purely from real-radiation corrections. Above the photon transverse- 
momentum cut of p\ > 100 GeV, the difference between the two distributions is very small. Below 
a missing transverse momentum of 100 GeV, the cross section is only non-vanishing starting from 
the NLO. Figure [ 6 ] shows a perturbative instability around pr,miss ~ 100 GeV. This instability 
originates from the incomplete cancellation of virtual and real corrections close to the phase space 
boundary (see Ref. [84] for a discussion of this phenomenon.) This class of singularities is inte- 
grable and does not alter the inclusive cross section, but would require a resummed computation 
to achieve a reliable differential prediction close to the boundary. 
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y/s = 7 Te V 


yfs = 8 TeV 


Neutrinos 

pif > 90 GeV 

pf > 100 GeV 

Photon 

p^ > 100 GeV 

p^ > 130 GeV 


|r/ 7 | < 2.37 


Frixione isolation with e 1 = 

o 

Oi 

o 

Jets 

pf > 30 GeV 


< 4.4 

rf et <4.5 


./Vjet > 0 (inclusive) or Aj et = 0 (exclusive) 

Separation 

AA( 7 ,jet) > 0.3 


Table 4: Event selection criteria for pp —» 


^ [TeV] 


OLO [fb] 

onlo [fb] 

Onnlo [fb] 

^ATLAS [fb] 

7 

A jet > 0 

Ajet = 0 

o9^-7+0.3% 

zu.z _Q # 9 % 

41 23 +4 ' 1% 

29 3fi +L2% 
zy.oo _ 1 3% 

46.01^5 

OQ Qtt + 1-0% 

28.85_ 09% 

±4.3 (stat) 
44.3 ±6.7 (syst) 

±1.7 (lumi) 

±3.3 (stat) 
38.7 ±4.3 (syst) 

±1.3 (lumi) 

8 

A jet > 0 
Ajet = 0 

14 n +1 - 1% 

Zt).DD_3_ 9% 

1 5 09 +1,6% 
io.uy _ 1 9% 

26.947“ 

14.893;“ 



Table 5: pp —> cross sections at LO, NLO and NNLO. Scale uncertainties are evaluated as 

in Table [21 The numerical uncertainty of the NNLO prediction from statistical error and finite 
r cut is conservatively estimated to be about 0.5%. The last column provides the corresponding 
result by ATLAS. 
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Figure 6 : Photon transverse momentum (left) and missing transverse momentum (right) distri¬ 
bution for pp —> vpvti a f LO (blue, dotted), NLO (red, dashed) and NNLO (green, solid). The 
lower panel shows the NNLO/NLO ratio. 


We can also study the transverse-mass distribution of the uv^ system, defined as 


(m 


VV') 

T 


2 



(5) 


Figure dsliows the transverse-mass distribution in the inclusive (left) and exclusive case (right). 
Transverse masses below m ^ 7 ~ 200 GeV are not allowed in LO kinematics and thus are only 
populated by real corrections starting from the NLO. This leads to an increased impact of the 
NNLO corrections in the region < 200 GeV in the inclusive case, with corrections of about 
100% compared to the NLO prediction. When applying a jet veto, this effect vanishes almost 
completey, indicating that relatively hard QCD radiation is necessary to overcome the LO kine¬ 
matics phase space constraint. In fact, with a 30 GeV jet veto present, the real radiation does 
only populate the phase space down to m"’ ! ~ 187 GeV at NLO and down to m ^ 7 ~ 175 GeV at 
NNLO. 

Figure [7] also shows the contribution from gluon fusion separately, which again is quite small 
and amounts to less than 2% of the fiducial cross section in the inclusive case and about 3% in 
the exclusive case. 


3.5 pp —> tvf) 

We now present results for pp —> at y/s = 7TeV and 8 TeV. We again use the event selection 

criteria adopted in the ATLAS analysis |9}. This set of cuts is identical to that used in the 
pp —> £ + £~^f analysis, apart from the fact that the cut on the invariant mass of the leptons is 
replaced with a cut on the missing transverse momentum (which coincides with the transverse 
momentum of the neutrino from the W decay) of > 35 GeV. As in the case of Z 7 in our 
y/s = 8 TeV analysis the rapidity of the jets is required to be |pj et | < 4.5. A summary of all cuts 
and event selection criteria can be found in Table EO 
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Figure 7: Transverse-mass distribution of the Z 7/77 system in the inclusive (left) and exclusive case 
(right) at LO (blue, dotted), NLO (red, dashed) and NNLO (green, solid). The loop-induced gluon 
fusion contribution is also shown (pink, dash-dotted). The lower panel shows the NNLO/NLO 
ratio. 


a/s = 7 Te V a/s = 8 TeV 


Lepton 

p e T > 25 GeV 
\rj\ < 2.47 

Neutrino 

p u T > 35 GeV 

Photon 

Pp > 15 GeV (soft pp cut) or p^ > 40 GeV (hard p^. cut) 

p 7 < 2.37 

Frixione isolation with e 1 = 0.5, R = 0.4, n — 1 

Jets 

anti-7'T algorithm with D = 0.4 
p J T et > 30 GeV 

|/ et | < 4.4 |/ et | < 4.5 

IVjet > 0 (inclusive) or A(j et = 0 (exclusive) 

Separation 

AR(£, 7 ) > 0.7 

Ai?(£/ 7 ,jet) > 0.3 


Table 6 : W ± (—y v^i) 7 cuts and event-selection criteria. 
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^ [TeV] 


c^LO [pb] 

ONLO [pb] 

CTNNLO [pb] 

catlas [pb] 

7 

Aj e t A 0 

Ajet = 0 

0.8726l®;®g 

2.05813 

1 7Q7 +5 - 2% 
i.oyj_ 58% 

2.45313 

1-49313 

±0.03 (stat) 
2.77 ±0.33 (syst) 
±0.14 (lumi) 

±0.03 (stat) 

1.76 ±0.21 (syst) 

±0.08 (lumi) 

8 

Ajet > 0 
Ajet = 0 

0-989313 

2 4f)l +7 - 4% 
z.4ui_ 74% 

1 587 +5 ' 5% 

i.°o ( 3% 

2-8847“ 

1.6917“ 



Table 7: W ± (—> ug£) 7 cross sections with the soft p^ cut (p\ > 15GeV). Scale uncertainties 
are computed as in Table [2j The numerical uncertainty of the NNLO prediction from statistical 
error and finite r cut is conservatively estimated to be about 0.8%. The last column provides the 
measured cross sections provided by ATLAS. 


yfi [TeV] 


OLO [lb] 

0NLO [lb] 

<7NNLO [fb] 

7 

CD 

IV 

0 

117 C+2.6% 

iio.o_ 37% 

707 q+9-0% 

oao.y_ 73% 

497 1 +5 - 3% 

±yi.i_ 47% 

8 

a> 

IV 

0 

1 n+ 3 . 5 % 

100.U4 6% 

478.613 

fi04 q+ 5 - 2% 
ou±.o_ 45% 


Table 8 : W ± (—> U{;£) 7 cross sections with the hard p 1 cut > 40 GeV). Scale uncertainties are 
computed as in Table [2j The numerical uncertainty of the NNLO prediction from statistical error 
and finite r cut is conservatively estimated to be about 0.5%. 


All results in the following will be presented summed over the W charges, i.e. we combine 
the processes pp —>■ W + ^ and pp —>■ W~ 7 , to facilitate the comparison with experimental data. 
The predicted fiducial cross sections both for the inclusive and the exclusive case can be found 
in Table [71 In the inclusive case, the NLO corrections are quite large, and amount to about 136- 
143%. The NNLO corrections increase the NLO result by 19-20%. The impact of higher order 
corrections is thus much larger than in the case of Z 7 production. We will come back to this point 
in Section [3761 

Table [7] also shows the cross sections measured by ATLAS. The measurement of the inclusive 
cross sections shows a 2 a excess with respect to the NLO prediction, which is reduced to well 
below 1 a when including the NNLO corrections. 

The impact of QCD corrections at NLO and NNLO is reduced to 60% and 7%, respectively, 
when the jet veto is applied (Aj et = 0). As discussed in Section 1331 such an effect is expected 
and apparently leads to a more stable perturbative prediction, but also to the possible need of 
more conservative procedures to estimate perturbative uncertainties. In the exclusive case, the 
excess of the measured fiducial cross sections over the theoretical prediction is reduced from 1 . 6 cr 
to 1.2 a when going from NLO to NNLO. We note that the scale variations at NLO significantly 
underestimate the impact of the NNLO corrections, in particular in the inclusive case. 

Figure [HI shows the photon transverse-momentum distribution in comparison with the ATLAS 
measurement, both in the inclusive and in the exclusive case. Although the experimental uncer- 
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Figure 8 : Photon transverse momentum distribution in the inclusive (left) and exclusive case 
(right) at NLO (red, dashed) and NNLO (green, solid) compared to ATLAS data. The lower panel 
shows the data/theory ratio. In the upper panel, only experimental uncertainties are shown. The 
lower panel also shows theoretical uncertainty estimates from scale variations. 


tainties are large, the agreement between data and theory is clearly improved when including the 
NNLO corrections, in particular if no veto on jets is applied. 


Figure E] shows the W 7 cross section differential in the transverse mass of the £277 system, 
normalized by the total fiducial cross section at the respective order. The transverse mass is 
defined here as 






( 6 ) 


The calculation is done with the hard photon transverse-momentum cut p\ > 40 GeV. The 
corresponding fiducial cross sections at LO, NLO, and NNLO are reported in Table [H] The impact 
of QCD radiative corrections is 242-260% and 26% at NLO and NNLO, respectively. I 11 Figure [9] 
due to the normalization, the large overall corrections mostly cancel out, in particular at high 
transverse mass, and we observe only a slightly improved agreement with data when going from 
NLO to NNLO. 

The increased relative impact of NLO and NNLO corrections when a harder p^ cut ( p ^ > 
40 GeV) is applied can, in analogy to the Z 7 case (see Section lT3l) . be better understood by 
studying the transverse-mass distributions with the soft and hard pl^ cut in more detail. The 
corresponding plots with a finer binning are shown in Figure [TOl When p vp > 15 GeV, for Born 
kinematics the transverse mass has a lower bound at about m ^ 7 > 75 GeV, i.e. below the W —>■ 
£277 peak. When the photon transverse-momentum cut is increased to 40 GeV, the lower bound 
increases to rn^ 1 > 100 GeV, and the W —> £277 peak is only populated by real emissions starting 
from the NLO, leading to large corrections in the region where the cross section is sizeable, and 
thus explaining the effect on the fiducial cross section. 
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Figure 9: Transverse-mass distribution of the system, normalized to the respective fiducial 

cross section at NLO (red, dashed) and NNLO (green, solid), compared to ATLAS data. The 
lower panel shows the data/theory ratio. Only experimental uncertainties are shown. 



Figure 10: Transverse-mass distribution of the 7 system at LO (blue, dotted), NLO (red, 
dashed) and NNLO (green, solid) for p\. > 15 GeV (left) and > 40 GeV (right), in the inclusive 
case (Aj et > 0). The lower panel shows the NNLO/NLO ratio. 
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process 

PT,cut 

[TeV] 


^nlo/^lo 

cnnlo/ Unix) 

Z (-> £+£-)7 


7 

CD 

IV 

0 

+50% 

+ 8 % 


soft 


-^jet — 0 

+27% 

+3% 


8 

-^jet ^ 0 

+50% 

+ 8 % 




-^Yjet — 0 

+25% 

+3% 


hard 

7 

Njet ^ 0 

+79% 

+17% 



8 

-^jet ^ 0 

+82% 

+18% 

Z (-> vpi) 7 


7 

-^jet ^ 0 

+57% 

+ 12 % 



Njet = 0 

+ 12 % 

- 2 % 



8 

-Njet + 0 

+ 68 % 

+14% 




-^jet = 0 

+7% 

- 1 % 

W (-► ivg) 7 


7 

O 

Al 

a; 

+136% 

+19% 


soft 

{ 

Njet — 0 

+60% 

+7% 


a 

a> 

IV 

0 

+143% 

+ 20 % 




-Njet — 0 

+60% 

+7% 


hard 

7 

-^Yjet ^ 0 

+242% 

+26% 



8 

Njet > 0 

+260% 

+26% 


Table 9: Summary of the relative NLO and NNLO corrections in the channels under investigation, 
Z (—* £ + £~) 7 , Z (—> V{V{) 7 , and W ± (—> i^vg) 7 . Numbers are reported for both soft and hard 
cuts. 
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3.6 The difference between W 7 and Z 7 

It is interesting to compare the relative size of the NLO and NNLO corrections to the Z 7 and 
W 7 processes we have considered. The results are summarized in Table |9j It is clear that the 
IT 7 process features much larger radiative effects with respect to the Z 7 processes. This should 
be contrasted to what happens in the case of inclusive W and Z boson production, where QCD 
radiative corrections are essentially identical [85]. It is thus the emission of the additional photon 
that breaks the similarity between the charged current and the neutral current processes. 

Restricting the analysis to NLO for the moment, the main source for the difference between 
Z 7 and fT 7 can be traced back to the gq and gq channels, which contribute a moderate, negative 
amount to the cross section in Z 7 production, but are large and positive for W ± 7 - It turns out that 
this effect is driven by resonant W 7 contributions to the cross section, i.e. by pp —> W{—¥ ivf )7 
topologies, and not by pp —> W —> £(—> £ 7)17 topologies, where the photon is emitted from the 
final-state lepton. These two contributions can only be separated in a gauge-invariant way if the 
W bosons are treated as on-shell particles, i.e. in a narrow-width approximation. By studying the 
LO contributions to the Z 7 and W 7 cross sections (see Figure [Hand Figure [2]) it turns out that in 
W 7 there is an additional Feynman diagram in which the photon is radiated off the W boson (see 
Figure l2dlh This additional diagram is responsible for a radiation zero [ 86 ], an exact zero present 
in the on-shell partonic W 7 tree-level amplitude at cos0* = 1/3, where 6 * is the scattering angle 
in the centre-of-mass frame. This radiation zero gets diluted by the convolution with the parton 
densities and by off-shell effects, but it is responsible for the suppression of the Born level W 7 
cross section with respect to Z 7 . Real radiation appearing at NLO breaks the radiation zero, and 
thus the relative impact of higher-order corrections is significantly increased. 

To quantitatively test this effect we consider the pp —> £vf] and pp —?> £ + £~7 processes studied 
in Section 13.31 and Section 13.51 with the same selection cuts. Contrary to what was done in the 
previous sections we disable the contributions from final state radiation and use the narrow-width 
approximation. I 11 Figure [IT] (left) we plot the distribution in the rapidity difference A ^ 7 between 
the charged lepton and the photon 1ST]. 

We see that the LO distribution shows a pronounced dip at central rapidities. Although diluted 
by the convolution with the parton densities, the dip is clearly visible, and is responsible for the 
suppression of the W 7 cross section. Since real radiation does not respect the radiation zero, the 
dip is filled up by radiative corrections. Roughly speaking, the NLO is a de facto LO prediction 
in the region of the dip and the NNLO corrections are thus relatively large as well. In contrast to 
IF 7 , the Z 7 amplitude does not exhibit a radiation zero and, consequently no dip appears in the 
rapidity-difference distribution, as can be seen in Figure [Til (right). 

The presence of the radiation zero, and of the corresponding dip in the A ^ 7 distribution, are 
thus the reason for the increased importance of radiative corrections to the W 7 processes. 


4 Summary and discussion 


In this paper we have presented the first complete and fully differential computation of QCD 
radiative corrections to W 7 and Z 7 production at hadron colliders. More precisely, we have 
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Figure 11: Rapidity difference between the charged lepton and the photon for W 7 (left) and Z 7 
production (right) at LO (blue, dotted), NLO (red, dashed) and NNLO (green, solid). The lower 
panel shows the NNLO/NLO ratio. Final-state radiation has been disabled for these plots. 


considered the processes pp —$■ £ + £~ 7 , pp —$■ and pp —> iu£j, where, in the first case, 

the lepton pair £ + £~ is produced either by a Z boson or a virtual photon. The diagrams in 
which the photon is radiated off the final-state charged leptons were consistently included. We 
have presented quantitative predictions for fiducial cross sections for pp collisions at y/s = 7 
and 8 TeV, and for various kinematical distributions (only at -y/s = 7 TeV). The impact of QCD 
radiative corrections strongly depends on the applied cuts. In the case of Z 7 , the impact of NNLO 
corrections is generally moderate, ranging from 8 % to 18%. We have also shown that the loop 
induced gluon fusion contribution is generally small, and it accounts for less than 10 % of the full 
0(a$) correction. In the case of W 7 production the NNLO effects are more important, and range 
from 19% to 26%. The larger impact of QCD radiative effects in the case of W'j production is a 
well known consequence of a radiation zero [ 86 ] existing in the W 7 amplitude at Born level. This 
effect produces a suppression of the LO distribution in the rapidity difference between the charged 
lepton and the photon, and NLO and NNLO corrections are thus quite significant. As expected, 
the impact of QCD radiative effects is strongly reduced when a jet veto is applied (Aj et = 0), 
being smaller than 3% in the case of Z 7 , and about 7% in the case of IR 7 . 

We add few comments on the remaining uncertainties affecting our NNLO results. The uncer¬ 
tainties from missing higher-order contributions were estimated through scale variations, which are 
performed through independent variations of the renormalization and factorization scales around 
their central value (without constraints on their ratio). In the inclusive case the NNLO scale 
uncertainties obtained in this way are of the order of ±(1 — 2 )% in the case of pp —>■ £ + £~7 (see 
Table [2]), ±(2 — 3)% in the case of pp —$■ Z 7 I 77 (see Table [5]), and ±4% in the case of pp —>■ £vff 
(see Table [7D . The comparison of the NNLO predictions to what is obtained at NLO shows that 
the NNLO-NLO difference is larger than the NLO scale dependence. We thus conclude that, as 
usual, scale variations can give only a lower limit on the true perturbative uncertainty. However, 
the NNLO is the first order at which all partonic channels are accounted for, and we believe that 
the NNLO scale uncertainties obtained in the case AQt > 0 should provide the correct order of 
magnitude of the true uncertainty. As discussed, the situation is different for the case Aj et = 0 , 
in which the scale uncertainties are even smaller. Most likely, in this case a more conservative 
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approach has to be adopted to obtain a realistic estimate of the perturbative uncertainty (see 

e.g. mm)- 


The other source of uncertainty affecting our perturbative QCD calculations is the one coming 
from the PDFs. The PDF uncertainties at 68 % CL that we obtain on our fiducial cross sections 
are at the 1% — 2% level, both at NLO and NNLO, and are thus of the same order, or smaller, 
than the perturbative uncertainties. We have checked that, by using CT10 ( 88 ] and NNPDFs [89] . 
the differences we obtain with the default MMHT result are of the same order. 

The quantitative predictions we have presented for i/s = 7 TeV were obtained by using the 
same cuts adopted by the ATLAS collaboration in their measurement of the hFy and Z 7 cross 
sections [91]. We have presented a comparison to ATLAS data, both for the fiducial cross sections 
and for some kinematical distributions. In the case of Z 7 production the NNLO corrections 
slightly improve the agreement with the data, which, however, have large uncertainties. The only 
exception is the case pp —> with a jet veto (Aj et = 0), for which a 1.7 a discrepancy with the 
ATLAS result remains. In the case of W 7 , the ATLAS result overshoots the NLO prediction by 
about 2a. The large NNLO corrections reduce this excess to below la. The NNLO corrections 
improve the agreement with the data also for the kinematical distributions we have studied, in 
particular for the pl^ distribution. However, it is known that this distribution, at large p^, is 
affected by sizeable effects from EW corrections [23H25] . The impact of EW corrections depends 
on the way the additional photon radiation is treated. It is thus difficult at present to draw 
conclusions on the data/theory agreement in the high-p^ region. More generally, we think it will 
be important to consistently combine the QCD calculations presented in this paper with a NLO 
computation of EW correction, which is, however, left for future work. 
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